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Abstract 

This document reviews the definition of the kernel distance, providing a gentle introduction tailored to a 
reader with background in theoretical computer science, but limited exposure to technology more common 
to machine learning, functional analysis and geometric measure theory. The key aspect of the kernel distance 
developed here is its interpretation as an L 2 distance between probability measures or various shapes (e.g. 
point sets, curves, surfaces) embedded in a vector space (specifically an RKHS). This structure enables several 
elegant and efficient solutions to data analysis problems. We conclude with a glimpse into the mathematical 
underpinnings of this measure, highlighting its recent independent evolution in two separate fields. 
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G\ • 1 Definitions 

Let K : R d x R d — > M be a similarity function with the property that for any x,K(x, x) — 1, and as the distance 
between x and y increases, K(x, y) decreases. A simple example of a kernel is the Gaussian kernel K(x, y) = 

o ; 

Definition 1.1 (Kernel Distance [3,4,6-8, 15, 16]). Let P and Q be sets of points in R d . The kernel distance 
^ | between P and Q is 



P^Pp'eP qeQq'eQ pePqeQ 



Note that the kernel distance is defined in terms of its square. We will continue this formulation so as 
to avoid the repeated use of square roots. While D K (P,Q) satisfies symmetry and one side of the identity of 
indiscernables, it might not in general be a metric, or even a pseudometric. In Section 4 we will discuss specific 
conditions on the similarity function K for which D K is a metric; in short, K must be a positive definite kernel. 

> : 
X 

H ■ 2 From Similarities to Distances 

C3 ■ 

The construction of the kernel distance involves a transformation from similarities to distances. This transfor- 
mation takes the following general form. Given two "objects" A and B, and a measure of similarity between 
them given by K(A,B), then the induced distance between A and B can be defined as the difference between 
the self-similarities K(A,A) + K(B,B) and the cross-similarity K(A,B): 

d{A,B) = K{A,A) + K{B,B)-2K{A,B). 

(Note that the multiplicative factor 2 is needed to ensure that d{A,A) = 0.) This construction can also be 
motivated set- theoretically: we can express the cardinality of the symmetric difference SAS' between two sets 
S,S' (a measure of the distance between them) as the expression |SAS'| = \S\ + \S'\ — 2\S C\S'\, where now the 
cardinality of the intersection plays the role of the similarity function K. 



*This is a "gentle" introduction, which means that where necessary, it will sacrifice mathematical rigor for ease of exposition. The 
cited references provide a more rigorous treatment of the material. 
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2.1 Point-wise Similarities 



Consider two points p,q. The kernel distance D^({p},{q}) = K(p,p) +K(q,q) - 2K(p,q) — 2(1 - K(p,q)) 
(when K{p,p) = 1). Here, the expression 1 —K{p,q) acts like a (squared) distance between the two points. 
Alternately, we can view the kernel distance as given by the sum of the self-similarities of p and q, from which 
we subtract the cross-similarity K(p, q). Note that 

• if the two points are identical, this expression is zero, i.e x = y => D K (x,y) — 0. 

• D K (p,q) = D K (q,p). 

• As the points become more distant, the self-similarities remain the same while the cross-similarity de- 
creases, which increases D K . 

Formulating the kernel distance in terms of self- and cross-similarities allows us to recover the general form 
of the definition. Let the cross-similarity of point sets P and Q be given by k(P,Q) = Xi pe p 2q 6 Q^(P> <?)• Tnen 
we can write 

D 2 K {P, Q) - k(P, P) + k(Q, Q) - 2k(P, Q). 
2.2 Comparing data with uncertainty 

Another way of interpreting the kernel distance comes from modeling data with uncertainty. If we have two 
sets of points P and Q, a naive way of estimating their difference would be to compute the symmetric difference 
|PAQ| = |P| + |Q| — 2|PnQ|. Of course, in all but a few cases, this number would be meaningless, since |P nQ| 
would equal zero unless points exactly coincided. We can replace the sharp similarity function k(p, q) = l p=q by 
the smoother similarity function K(p, q), where K(p, •) represents uncertainty in the location of p (for example, 
as a Gaussian). In this case, the above symmetric difference \PAQ\ then turns into D\(p,q). 



3 Generalizations 

One of the features of the kernel distance is that it can be generalized beyond point sets to distributions in 
space and even to higher dimensional geometric structures like curves and surfaces, where it is often called 
the current distance [16]. In this section we merely define the appropriate generalizations, giving some simple 
intuition where appropriate. A deeper treatment of the mathematics underpinning these constructions will be 
presented in Section 5. 



3.1 Measures and Distributions 

We can think of a point set as a set of 5-functions at each location, which gives us a representation of the point 
set as a measure over R d . We can then weight these functions, which allows us to compare weighted point sets. 
Let 7 = (P, w) be a point set with an associated weight function w : P — * R, and let Q = (Q, w') be another 
such weighted point set. We can define the cross-similarity k(!P, Q) — T^^^-p ^~!g^-r) w(p')K(p, q^)w' '(<?) and then 
as before define 

D*CP, Q) = k(9, IP) + k(Q, Q) - 2k(9, Q). 

In general, given two distributions over R d , we can define the kernel distance between them merely by replacing 
the sum in the cross similarity by an integration. Specifically, if we are given two distributions /i, v, we define 
the cross-similarity as 



after which D K can be defined in the usual manner. 



K(p,q)dn(p)dv(q) 
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3.2 Curves 



Let us now consider the case when P and Q are curves. For convenience we will assume that they are C 1 
continuous (i.e. continuous curves with well-defined tangent vectors everywhere). While it is possible to define 
the kernel distance between them by treating the curves as (infinite) sets of points, this approach ignores the 
fact that we would like to match not only the locations of points but the curve gradients. 

Let tp(p) denote a unit tangent vector at the point p e P on the curve. We can then define the pointwise 
similarity between the two points p e P and q e Q as the function k(p,q) = K(p, q){t P {p), tq(q)). The cross- 
similarity between P and Q can then be defined as 



<P,Q) = 



K(p,q)(t P (p), t Q (q)>. 



(3.1) 



JP JQ 



3.3 Surfaces 

The same construction as above can be applied to orientable surfaces of co-dimension 1 in R d , with the differ- 
ence being that we use a unit normal at a point (which is a vector), in place of its tangent. 



4 Reproducing Kernel Hilbert Spaces 

Not all similarity functions K yield a distance metric. Consider the box similarity given by 



K(x,y) 



1 ||x-y||<2 
otherwise. 



Let A be a set of three points all lying close to the origin in R 1 . Fix B to be the set {(—1 — e), (1 + e), (0)}. It is 
easy to verify that D^(A,B) — —2, and thus D K (A,B) is not even defined. 
There is a simple sufficient condition on K to ensure that D K to be a metric: 

K must be a positive definite kernel. 

For those well-versed in the theory of reproducing kernel Hilbert spaces [1], this observation follows easily by 
replacing the K(x,y) terms by the appropriate inner products in the feature space. For those unfamiliar with 
this theory, an explanation follows. The material here is cobbled together from Daume's[2] and Jordan's [9] 
notes on RKHSs. 



4.1 Positive Definite Matrices 

Positive definite kernels generalizes the idea of a positive definite matrix. Let us start with a simple matrix 
kernel: K{x,y) = x J Ay. If we set the matrix A to be identity, we get the standard Euclidean inner product 
(x,y). In general however, we can set A to be any positive definite matrix and guarantee that the resulting 
function x T Ay satisfies the properties of an inner product. 



A lifting map. Any positive definite matrix A can be decomposed as A = B J B. Substituting this back into the 
bilinear form x T Ay, we get the expression x T B T By, which we can rewrite as (Bx,By). If we now define the 
linear map $(x) = Bx, we can rewrite the original matrix kernel as 

There is further structure in the mapping $(x) that we can explore. Since K is symmetric as well as real-valued, 
it can be written in the form K — QAQ T where the columns of Q are eigenvectors of K, and A is a diagonal 
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matrix consisting of the eigenvalues of K. Since K is positive definite, all the elements of A are positive, and 
we can thus set B = A 1/,2 Q T , which means that each coordinate of <3>(x) is of the form -JH{v, x), where v is an 
eigenvector of K and A is the associated eigenvalue. More generally, this means that the eigenvectors of K form 
an orthonormal basis describing the vector space $(x). 

4.2 General Positive Definite Kernels 

If K is a general function (i.e not a matrix), much of the above theory carries over in a similar manner, where 
instead of vector spaces we now reason about function spaces. 

A note on function spaces. 

While not entirely rigorous, it is helpful to think of function spaces as merely infinite-dimensional vector spaces. A 
function / : X — » R can be thought of as a vector whose "dimensions" are indexed by elements of X, so that /(x) = c can 
be translated as "the xth coordinate of the vector / is c." Similarly, the function K : X x X — > R can be interpreted as a 
matrix indexed by elements of X where the (x, y)th entry is K(x,y). 

Integration can then be viewed as a product. For example, the integral J K(x,y)f(y)dy can be interpreted as a general- 
ized matrix product between the "matrix" K and the "vector" / . We will also look at "eigenfunctions" as the generalization 
of eigenvectors. 

We first define a positive definite kernel. 

Definition 4.1. A function K : X xX — > R is a positive definite kernel if for any n and any set {x 1; x 2 , . . . , x n } c X, 
the matrix A — (a i; - = K(x u x ; )) is positive definite. 

Any positive definite kernel induces a Hilbert space (the reproducing kernel Hilbert space) 'K. There are two 
different constructions yielding JC. We will focus here on the method based on Mercer's theorem, since that 
yields an orthonormal representation of that will be useful. 

A matrix K defines the linear operator T K {f) — Kf , here operating on a vector / and producing another 
vector. Similarly, a positive definite kernel K defines a linear operator on functions given by 

[T K ifm = jK(;y)f(y)dy. 

Again in this continuous setting, T K is an operator that now takes a function / as input and returns a function 
[T K (J )](•)• The operator T K is linear by virtue of the linearity of the integral, and so it has eigenvalues and 
eigenfunctions. 

Mercer's theorem can then be seen as the analog of the above decomposition of the positive definite matrix 
A as B J B,B = A 1 ^ 2 Q T . 

Theorem 4.1 (Mercer). IfK is a continuous symmetric positive definite kernel, then there is an orthonormal basis 
|v ; } consisting of eigenfunctions ofT K (with associated eigenvalues k { ) such that 

00 

K(x, y)-^ A ; Vj(x)Vj(y). 
o 

The trick here is in thinking of V;(x) as the xth "coordinate" of the "eigenvector" v t . 

The main consequence of this result is that we can now describe "K as the vector space spanned by the v h 
with the associated inner product given by 

By scaling the coordinates c ; by we can use the Euclidean inner product instead, giving us a Euclidean 
space. We summarize with the following statement, which captures the key properties of positive definite 
kernels. 
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Theorem 4.2. For any positive definite kernel K : X x X — * R, there exists a Euclidean space 'K and a lifting map 
$ : X -> IK such that 

K(x,y) = <<K*),<Ky)>. 
4.3 The Kernel Distance as a Hilbertian metric 

We can now rewrite the kernel distance D K in terms of the associated lifting map Using the linearity of the 
inner product, it is easy to show that 

D K (P,Q) = \\Y i Hp)-Y i ^ ■ 

pep qeQ 

This of course immediately implies that D K is a pseudometric. A further technical condition on K is needed to 
ensure that D K (P,Q) = => P — Q; the reader is referred to [14] for more details. 
There are four important consequences of this recasting of the kernel distance. 

1. The kernel distance embeds isometrically in a Euclidean space. While in general "K might be infinite 
dimensional, the Hilbert space structure implies that for any finite collection of inputs, the effective 
dimensionality of the space can be reduced via projections to a much smaller finite number [10, 12, 17]. 

2. Most analysis problems are "easier" in Euclidean spaces. This includes problems like near-neighbor find- 
ing and clustering. The embedding of the kernel distance in such a space means that we now have access 
to a number of tools for analyzing collections of shapes. Indeed, a complicated shape can be represented 
as a single vector f>(P) = Y* p eP *0>) in the RKHS [13]. 

3. The embedding "linearizes" the metric by mapping the input space to a vector space. This is indeed the 
primary reason for the popularity of kernel methods in machine learning. In the context of shape analysis, 
it means that many problems in the analysis of shape (finding consensus, averages, etc) can be solved 
easily by exploiting the linear structure of the lifted space. 

4. The complexity of computing the kernel distance is reduced significantly. If we assume that "K is approx- 
imated to within the desired error by a space of fixed dimension p, then in this space, computing the 
kernel distance between two point sets of total size n takes O(np) time, instead of G(n 2 ) time. Since p 
will in general be logarithmic in n (or even independent of n), this is a significant improvement. 

5 Mathematical Motivation 

There are two distinct motivations that yield the kernel distance as formulated in (1.1). The first motivation 
is based on the desire to metrize distributions, to construct a metric on distributions such that a convergent 
sequence with respect to this distance metric also converges in distribution. The second motivation comes from 
shape analysis in an attempt to use a geometric measure-theoretic view of shapes to define a distance between 
them that does not rely on explicit correspondences between features of shapes. 

Both of these approaches rely on the idea of "test functions" and "dual spaces", which we explore next. 

5.1 Dual Vector Spaces And Test Functions 

Let V be a vector space, and let / be a continuous linear functional over V. In other words, / is a linear 
function that takes an element v e V, and returns a scalar /(v). The space of all such / is itself a vector space, 
and is called the dual of V, denoted V*. A very simple discrete example has V and V* represent the columns 
and rows, respectively, of a fixed matrix. We can use structures defined on V* to construct similar structures on 
V; this is particularly useful when V is an ill-formed space, but V* has more structure. 
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Consider another simple example. Let V be R d , in which case the dual space V* is also R d . Suppose we wish 
to define a norm on V. We can do this via pullback from a norm on V* by the following construction: 

||v|| = sup |w(v)|. 
weV*,||w||<l 

It is not hard to verify that this indeed satisfies the properties of a norm. In fact, if we set the norm on V* to be 
an (. p norm with p > 1, then the resulting norm || • || on V is £ q , where l/p + l/q = l (the dual norm). 

Notice that choosing different subsets of V* (by changing the norm) allows us to generate different norms 
on V. The case of i 2 is particularly instructive. Here, the set of all w such that ||w|| < 1 is the set of all vectors 
in the unit ball, and since |w(v)| is increased by making w larger, we can assume without loss of generality 
that ||w|| = 1, which means that we are considering all directions in R d . In this case, w(v) is merely (w, v), the 
projection of v onto the direction w, and so the maximization over w returns precisely the l 2 norm of v. 

Norms yield metrics in the usual way, by setting d(x,y) = ||x — y\\. Continuing the same construction, we 
can define a metric on V by the construction 

d(v,v') = sup |w(v) - w(v')|. 
wev*,||w||<i 

Again returning to our example in l 2 , this has a nice geometric interpretation. Since each w is a direction in 
R d , the expression being maximized is merely the length of the vector vv' when projected onto w. Clearly, 
maximizing this yields the actual l 2 distance between v and v' . 



5.2 Metrizing Measures 



With the above construction in mind, we can now build a metric over probability measures. The general 
construction, deemed an integral probability measure by Muller [11], takes two probability measures P and Q 
over a space M, and defines the distance 1 



fdP- fdQ 



dgr(P,Q) = SUp 

feS 

Here, the class 3~ is the dual space to the space of probability measures, and consists of a class of real-valued 
bounded measurable functions over M. We can think of J fdP as the action /(P) of / on P; this is linear since 
the integral operator is linear. 

Many interesting distances between distributions can be generated by picking a particular subset of 3~. Some 
well-known cases: 



If?={/| 



< 1}, then dj(P,Q) is the l x distance between P and Q. 



L =sup{ f -^f,x^y}, then the 



• If M, the domain of /, is endowed with a metric p, and we define 

resulting distance is the Kantorovich metric, also known as the Earth mover's distance [5]. 

Finally, and this is the case of relevance here, we can set 3 = {f \ ||/||jt - 1}> where "K is a reproducing 
kernel Hilbert space. In this setting, it can be shown that the resulting distance can be rewritten as 



D K (P,Q) = 



K(-,x)dQ(x) 



•K 



K{-,x)dP{x)- 

which, after some expansion and simplification, yields the familiar form 

Df(P,Q)= K{x,y)dP{x)dP{y)+ K(x,y)dQ{x)dQ(y)-2 I I K(x,y)dP{x)dQ{y) 



1 Note that j fdP is a sparser notation for the more explicit representation of the intregral J^/(x)P(x)dx or J/(x)dP(x). 
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which brings us back to (1.1). 

Note that by choosing the measures P and Q appropriately, we can generate the different variants of the 
kernel distance described in Section 3. 

5.3 Currents: Distributions with geometry 

We now move to the problem of comparing shapes, specifically point sets, curves and surfaces. The key view- 
point we will take is to treat the shape S as a probability measure. For point sets, this approach is easy to see, 
and for curves and surfaces, an appropriate discretization yields the desired probability measure. Once this is 
done, we can invoke the techniques from Section 5.2 to define a distance metric between two shapes. 

For point sets, this construction is sufficient; however, curves and surfaces have more information than merely 
the spatial location of points. The well known example comparing the Hausdorff distance to the Frechet dis- 
tance for comparing curves illustrates that a proper comparison of curves requires not only spatial information, 
but also information about orientation, and this continues to be true for surfaces. 

Geometric measure theory is the discipline that develops the tools to do this, and develops the idea of a cur- 
rent as the geometric measure-theoretic analog of a surface that elegantly encodes the orientation information 
of the surface. The idea of normal vectors generalizes to higher dimensions with k-forms. 

Wedge Products, k-vectors, and k-forms 

The cross product of two vectors in K 3 is a well-defined function whose main characteristics are that v x v = and 
vxw + wxv = 0. The wedge product (denoted by A) generalizes the cross product to arbitrary vectors v e R d . It has 
the property that v A v = (which implies that v Aw = — wAv"). A simple k-vector is denned as the wedge product 
of k vectors: w = v x A v 2 ■ ■ ■ A v k , and the space of all linear combinations of k-vectors is denoted A fc (K d ). Note that 
A 1 (K d ) = R d . k-vectors generalize the idea of a tangent vector, by capturing in one expression a k-dimensional subspace 
of the tangent space at a point. 

A differential form can now be understand as a linear functional over the k-vectors. Formally, a k-form co is an assignment 
of a linear functional co x at each point x on the manifold. This linear functional co x is a linear mapping from A^(R d ) to 
R, and the space of all such functionals is denoted A fc (IR d ). Since k-forms and k-vectors are dual by construction, we will 
use the expression co{v~) or (w,v) to denote the induced action of one on the other. 

"This can be shown by expanding (w + v) A (w + v) = 0. 

We can integrate a k-form co over a k-dimensional oriented manifold S. Using dS to denote the appropriate 
"tangent space displacement" (a generalized tangent vector) on the surface, we can write this integral as 
J" co(dS) or J" codS when the action is understood. By integrating on the k-vector dS at each point peS, 
instead of just the point p, we capture the orientation information of the manifold S, not just the spatial 
information. By duality, this can be seen as an action of S on co, and leads us (finally) to the definition of a 
current: 

Definition 5.1. A k-current is a continuous linear functional over the space of k-forms. 

Currents generalize oriented manifolds. They allow us to account for irregularities and sharp edges in a 
manifold, at the cost of including elements that may not look like standard manifolds. Very crudely, the 
relation between currents and oriented manifolds is akin to the relationship between (Schwarz) distributions 
and probability measures. 

We can now define a norm on the space of currents much as before, by placing a norm on the space of 
k-forms and pulling it back. The equivalent of the total variation norm (by computing supa>(v), WcoW^ < 1) 
is called the mass norm; unfortunately, much like the total variation, it does not capture variations between 
currents in a useful manner. For example, given two curves with disjoint support, the distance induced by the 
mass norm is 2, irrespective of the actual shape of the curves. This is similar to how the total variation distance 
between two distributions is always 2 if they have disjoint support. 
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By setting the space of test functions to be the unit ball in an RKHS as before, we retrieve the kernel distance 
as described in Section 3. For instance, consider two curves S and T so the -vectors dS and dT are the tangent 
vectors t s (-) and t T (-)- We can expand 



D K (S, T) 



K(-,x)dS(x)- K(-,x)dT(x) 



■K 



as 



D Z K {S,T)= (dSMYK(x,y)dS(y) + 



{dT{x)) l K{x,y)dT{y)-2 (dS(x)) J K(x , y)dT '(y) 



(t s (x)yK(x,y)t s (y)+ (t T (x)Y K(x,y)t T (y) - 2 (t s (x))^(x, y)t T {y), 
sJs JtJt JsJt 

where (-) T represent a vector transpose. Since K(x,y) is a scalar, we can factor it out of each term 



D*(S,T)=\ K(x,y)(t s (x),t s (y))+ \ K(x, y)(t T (x), t T (y)) - 2 \ K(x,y)(t s (x), t T (y)), 
JsJs JtJt JsJt 

precisely as in the similarity term from (3.1). This derivation of D K has led to the alternate name current 
distance. 
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